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Abstract 

We study the strong interactions among the heavy bottom spin-1/2 S;, baryon, 
nucleon and B meson as well as the heavy charmed spin-1/2 Sc baryon, nucleon and 
D meson in the context of QCD sum rules. We calculate the corresponding strong 
coupling form factors defining these vertices by using a three point correlation func¬ 
tion. We obtain the numerical values of the corresponding strong coupling constants 
via the most prominent structure entering the calculations. 
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1 Introduction 


In the recent years, snbstantial experimental improvements have been made on the spec¬ 
troscopic and decay properties of heavy hadrons, which were accompanied by theoretical 
studies on various properties of these hadrons. The mass spectrum of the baryons con¬ 
taining heavy quark has been studied using different methods. The necessity of a deeper 
understanding of heavy flavor physics requires a comprehensive study on the processes of 
these baryons such as their radiative, strong and weak decays. For some of related studies 
one can refer to references [1-10]. 

The investigation of the strong decays of heavy baryons can help us get valuable in¬ 
formation on the perturbative and non-perturbative natures of QCD. The strong coupling 
constants dehning such decays play important role in describing the strong interaction 
among the heavy baryons and other participated particles. Therefore, accurate determina¬ 
tion of these coupling constants enhance our understanding on the interactions as well as the 
nature and structure of the participated particles. The present work is an extension of our 
previous study on the coupling constants QAbNB and Qk^nd [Hj- Here, we study the strong 
interactions among the heavy bottom spin-1/2 S;, baryon, nucleon and B meson as well as 
the heavy charmed spin-1/2 baryon, nucleon and D meson in the context of QCD sum 
rules. In particular, we calculate the strong coupling constants and gj^^ND- These 

coupling constants together with the Qa^nb and QAcND discussed in our previous work, may 
also be used in the bottom and charmed mesons clouds description of the nucleon which 
can be used to explain the exotic events observed by different Collaborations. In addition, 
the determination of the properties of the B and D mesons in nuclear medium requires 
the consideration of their interactions with the nucleons from which the and are 
produced. Therefore, to determine the modihcations on the masses, decay constants and 
other parameters of the B and D mesons in nuclear medium, one needs to consider the 
contributions of the baryons Sb[c] together with the Afe[c] and have the values of the strong 
coupling constants g-s^NB and g-E^ND besides the couplings gA^NB and gA^ND [12-15]. In 
the literature, one can unfortunately hnd only a few works on the strong couplings of the 
heavy baryons with the nucleon and heavy mesons. One approximate prediction for the 
strong coupling gAcND was made at zero transferred momentum squared [4]. The strong 
couplings of the charmed baryons with the nucleon and D meson were also discussed in [7] 
in the framework of light cone QCD sum rules. 

This paper is organized in three sections as follows. In the next section, we present the 
details of the calculations of the strong coupling form factors among the particles under 
consideration. In section 3, the numerical analysis of the obtained sum rules and discussions 
about the results are presented. 

2 Theoretical framework 

This section is devoted to the details of the calculations of the strong coupling form factors 
and from which the strong coupling constants among the participating 

particles are obtained at = —q^ = —subsequently. In order to accomplish this 
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purpose, the following three-point correlation function is used: 

n = J d^y {0\T {vN{y) Vb[d]{ 0) f]^,[j:,]{x)) |0), (1) 

whith T being the time ordering operator and q = p — p' is the transferred momentnm. 
The currents tjn, Vb[d] and presented in Eq. (1) correspond to the the interpolating 

currents of the N, B[D] and respectively and their explicit expressions can be given 
in terms of the qnark held operators as 

hsqs,](x) = £pfc(^n*'^(x)C'7/,d-^(x))757^6[c]^(x), 

VN{y) = 

Vb[d]{0) = h(0)75b[c](0), (2) 

where C denotes the charge conjngation operator; and i, j and k are color indices. 

In the course of calculation of the three-point correlation function one follows two differ¬ 
ent ways. The hrst way is called as OPE side and the calculation is made in deep Euclidean 
region in terms of quark and gluon degrees of freedom using the operator product expan¬ 
sion. The second way is called as hadronic side and the hadronic degrees of freedoms are 
considered to perform this side of the calculation. The QCD sum rules for the coupling 
form factors are attained via the match of these two sides. The contributions of the higher 
states and continnum are suppressed by a double Borel transformation applied to both sides 
with respect to the variables p‘^ and 

2.1 OPE Side 

For the calculation of the OPE side of the correlation fnnction which is done in deep 
Euclidean region, where —)■ — cxd and —)■ —oo, one puts the interpolating currents 
given in Eq. (2) into the correlation fnnction, Eq. (1). Possible contractions of all quark 
pairs via Wick’s theorem leads to 

f I I 

X ^l5duS^J{y - x)'yf,CSi^^{y - x)C'7^S'“'^(p)755qc](-a^)7/.75 

- l5'luSa\y - x)'yi,CS^^^{y - , (3) 

where 5'6[c](x) and S'„[rf](x) are the heavy and light quark propagators whose explicit forms 
can be fonnd in Refs. [11,16]. 

After some straightforward calculations (for details refer to the Ref. [11]), the correlation 
function in OPE side comes out in terms of different Dirac structures as 

= ni(g2)^5 + n2(g2)^75 + n3(g')#^75 + n4(g2)#75. (4) 


2 



Each nj(g^) function involves the perturbative and non-perturbative parts and is written 
as 



= ds ds 




(s — p‘^){s' — p'^) 


(5) 


The spectral densities, pi{s, s', q"^), appearing in Eq. (5) are obtained from the imaginary 
parts of the 11* functions, i.e., pi{s,s',q‘^) = -Im[Ili\. Here to provide examples of the 
explicit forms of the spectral densities, among the Dirac structures presented above, we 
only present the results obtained for the Dirac structure *^ 75 , that is s',q'^) and 

^non-pert^^^ s', q^), which are obtained as 




pi pl—x 

dx dy 
'0 Jo 


327r‘^(x + y — ly 


— y — 2x + Sxy^ 


- 3ml[^-^x{x + y - 1) - mu{x + Ay - 2) + md{3x + 6y - 2) 


q m!,[c]X 


X 


y -Sy"^ + x^iSy - 1) + x{l -Qy + 8y^) + mqc](a; + y - 1) 
sx{8x^ — 3y — hx + 8xy) + s'(^8xy‘^ — hxy — 2y‘^ — x'^ + 8x^y^ 

— {x + y — 1) 3mdS (^Ax^ + y — y"^ — 7x‘^ + 12a:^^ + 3a; — 12xy + 8xy‘^'^ 

+ 3mdq^y {Ax — Ax^ + y — 8xy'^ + 3mds'y ^3 — 7a: + 4a;^ — lly + I2xy + 8y'^ 

— rriuS Ax^ + y{3 — 2y) + a:^(20?/ — 13) + a;(9 — 27y + 16?/^) 

+ rrius'yi^ — lAx + Ax^ — 26y + 20xy + 16?/^ j 

+ m**g^?/(lla; - + 2?/- 16a;?/- 3)l0[L2(s,s',g^)l, (6) 


and 


non—pert / i 2 \ 

P2 {s,s,q ) = 


{uu) 


247r2(g2 _ 


q 2 mb[c] (3q^ - - 2s' 

+ (cts —)-. — n/ o-^ + 


TT 1927r2(g2 _ mLi)2 


(^3mb[c]md - 3mumd + 3ml - q^ + s - s''^ 

Li{s,s',q‘^) 


mQpuu) 


b[c]i 
1 pl—x 


0 


- {uu)) I dx I dy ^^^^J ^ —^0 


167r2(g2 _ 77^2 

L2{s,s',q^) 


0 ao 


(7) 


where 0 [...] stands for the unit-step function and Li(s, s',g 2 ) and L 2 {s,s',q'^) are dehned 
as 


Li{s,s',q^) = s', 

/ f 2 \ 2 2 f 2 f f / \ 

L 2 {s,s ,q ) = —-t- sa; — sa: + s y + q xy — sxy — s xy — s y . (8) 
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2.2 Hadronic Side 


On the hadronic side, considering the quantum numbers of the interpolating helds one place 
the complete sets of intermediate Sfe[Sc], B[D] and N hadronic states into the correlation 
function. After carrying out the four-integrals, we get 

J.HAD ^ (0 hiV I N{p')){0 I r]B[D] I ^[A>](g))(Sfe[S,](p) | I 0) 

X {N{p')B[D]{q)\E,[i:,]{p)) + --- . (9) 


In the above equation, the contributions of the higher states and continuum are denoted by 
■ ■ ■ and the matrix elements are represented in terms of the hadronic parameters as follows: 


(OhwIW)) = 
(Sfe(p) I I 0) = 

{0 \ Pbid] \ B[D]iq)) = 
{N{p')B[D]{q) I S,[S,](p)) = 


\NUNip',s'), 

As(.[Sc]hsqSc] {P, ■s), 

.^B[D]fBlD] 

*- , -’ 

5'S(.AfB[ScArD]'hv(P^ 'S')f75'f^EqEc] (P) 'S)- 


( 10 ) 


Here and AsqEq are residues of the N and Sf,[Sc] baryons, respectively, /b[d] is the 
leptonic decay constant of B[D] meson and Pt.^nbi'ScND] is the strong coupling form factor 
among Sfe[Sc], N and B[D] particles. Using Eq. (10) in Eq. (9) and summing over the spins 
of the particles, we obtain 


jlHAD ^ B[D] _ AjvAsqEc]^E6WB[EeWE>] _ 

mqq + {p^ - “ "^v)(9^ - ^%[d]) 

X |(m7vmsqE,] - m,|^[E,])75 + ("^EqE,] - m^) ^75+ 4 ^75 - "^sqE,] # 75 } 

+ ■■■ . ( 11 ) 


To acquire the hnal form of the hadronic side of the correlation function we perform the 
double Borel transformation with respect to the initial and hnal momenta squared. 


= ,2 ^^B[D]fB[D] AjvAsqEq^E,VB[EeiVD] 

mqc] + ruu {q^ - m|pj) 

X |(m7vmsqEe] - "^lqEd)75 + ("^sqEj - tun) ^ 75 + 4 ^75 - ^e^e,] 4lb] 

+ ••• , ( 12 ) 

where and are Borel mass parameters. 

As it was already stated, the match of the hadronic and OPE sides of the correlation 
function in Borel scheme provides us with the QCD sum rules for the strong form factors. 
The consequence of that match for ;4j5 structure leads us to 


gj^bNBlT.cNDliq^) 


[nibic] + rnu){q^ - m|r^,) 

— g ivp gM'^ ---i— - 

^B[D]/Bp]"^E(,[Ec]'^^(^^^S6pc] ~ ''^sqEj) 
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X 


fSO 


ds 


ds'e e 


' (mb[^]+mu+md)‘2 


' (2m„+mjj)2 


pert/ / 2\ \ non—pert/ / 2 

(h (s,S,g)+P2 


(13) 


where Sq and Sq are continuum thresholds in Sb[Sc] and N channels, respectively. 


3 Numerical analysis 

Having obtained the QCD sum rules for the strong coupling form factors, in this section, 
we present the numerical analysis of our results and discuss the dependence of the strong 
coupling form factors under consideration on = —q^. To this aim, beside the input 
parameters given in table 1, one needs to determine the working intervals of four auxiliary 
parameters M^, sq and Sq. These parameters originate from the double Borel trans¬ 
formation and continuum subtraction. The determination of the working regions of them 
is made on the basis of that the results obtained for the strong coupling form factors be 
roughly independent of these helping parameters. 

The continuum thresholds sq and Sq are the parameters related to the beginning of the 
continuum in the initial and hnal channels. If the ground state masses are given by m and m' 
for the initial and hnal channels, respectively, to excite the particle to the hrst excited state 
having the same quantum numbers with them one needs to provide the energies — m 
and For the considered transitions, these quantities can be determined from well 

known excited states of the initial and hnal states [17] which are roughly in between 0.1 GeV 
and 0.3 GeV. From these intervals, the working regions for the continuum thresholds are 
determined as 34.9[6.5] GeV^ < so < 37.4[7.6] GeV^ and 1.04 GeV^ < -Sq < 1.99 GeV^ for 
the vertex B\^cND]. 

Here, we shall comment on the selection of the most prominent Dirac structure to 
determine the corresponding strong coupling form factors. In principle, one can choose 
any structure for determination of these strong coupling form factors. However, we should 
choose the most reliable one considering the following criteria: 

• the pole/continuum should be the largest, 

• the series of sum rule should demonstrate the best convergence, i.e. the perturbative 
part should have the largest contribution and the operator with highest dimension 
should have relatively small contribution. 

Our numerical calculations show that these conditions lead to choose the structure 4 7^75 as 
the most prominent structure. In the following, we will use this structure to numerically 
analyze the obtained sum rules. 

Now, we proceed to present the Borel windows considering the selected structure. The 
working windows for the Borel parameters and are determined considering again 
the pole dominance and convergence of the OPE. By requirement that the pole contribution 
exceeds the contributions of the higher states and continuum, and that the contribution 
of the perturbative part exceeds the non-perturbative contributions we hnd the windows 
10 [2] GeV^ < < 20 [6] GeV^ and 1 GeV^ < < 3 GeV^ for the Borel mass parameters 
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Parameters 

Values 

mb 

(4.18 ±0.03) GeV[17] 

rric 

(1.275 ±0.025) GeV[17] 

ma 

MeV[17] 

mu 

2.3tol MeV [17] 

ms 

(5279.26 ± 0.17) MeV [17] 

mo 

(1864.84 ±0.07) MeV [17] 

mN 

(938.272046 ± 0.000021) MeV [17] 


(5811.3 ± 1.9) MeV [17] 

ms. 

(2452.9 ± 0.4) MeV [17] 

fs 

(248 ± 23e.p ± 25vub) MeV [18] 

Id 

(205.8 ± 8.5 ± 2.5) MeV [19] 


0.0011 ±0.0005 GeV® [20] 

Asi, 

(0.062 ± 0.018) GeV® [5] 

Ae, 

(0.045 ± 0.015) GeV® [5] 

{uu){l GeV) = {dd){l GeV) 

-(0.24 ±0.01)® GeV® [21] 

(^) 

(0.012 ±0.004) GeV^ [22] 

mg(l GeV) 

(0.8±0.2) GeV^ [22] 


Table 1: Input parameters used in calculations. 


of the strong vertex TjbNB[TjcND]. For these intervals, our results show weak dependence 
on the Borel mass parameters (see hgures 1-2). 


10 12 14 16 18 20 1.0 1.5 2.0 2.5 3.0 




Figure 1: Left: = 0) as a function of the Borel mass at average values of 

continuum thresholds. Right: gi^iNsiQ'^ = 0) as a function of the Borel mass M' at 
average values of continuum thresholds. 

Subsequent to the determination of the auxiliary parameters, their working windows 
together with the other input parameters are used to ascertain the dependency of the 
strong coupling form factors on Q^. From our analysis we observe that the dependency of 
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Figure 2: The same as figure 1 but for nd{Q‘^ — 0 ). 


-30 -20 -10 


10 20 30 




Figure 3: Left: gs^NsiQ'^) as a function of at average values of the continuum thresholds 
and Borel mass parameters. Right: gEcND^Q'^) as a function of Q‘^ at average values of the 
continuum thresholds and Borel mass parameters. 


the strong coupling form factors on is well characterized by the following £t function: 


6'E6AfB[EcAfD](Q^) — Cl exp 


91 

C2 


+ C3. 


(14) 


The values of the parameters Ci, C 2 and C 3 for Ti^NB and Ti^ND can be seen in tables 2 and 
3, respectively. Considering the average values of the continuum thresholds and Borel mass 
parameters we demonstrate the variation of the strong coupling form factors with respect 
to for the QCD sum rules as well as the htting results in hgure 3. The hgure indicates 
the truncation of the QCD sum rules at some points at negative values of Q^. It can be 
seen from the hgure that there is a good consistency among the results obtained from the 
QCD sum rules and £t function up to these points. The £t function is used to determine 
the value of the strong coupling constant at and the results are presented in 

table 4. The presented errors in these results originate from the uncertainties of the input 
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parameters together with the uncertainties coming from the determination of the working 
regions of the auxiliary parameters. 


Parameters 

ci(GeV-^) 

C 2 (GeV 2 ) 

C 3 (GeV-i) 

Valnes 

0.72 ±0.20 

18.03 ±4.87 

9.60 ±2.88 


Table 2: Parameters appearing in the £t function of the coupling form factor for Ti^NB 
vertex. 


Parameters 

ci(GeV-^) 

C2(GeV2) 

C3(GeV-i) 

Valnes 

-0.006 ± 0.002 

-6.78 ± 1.89 

3.79 ± 1.02 


Table 3: Parameters appearing in the £t function of the coupling form factor for Ti^ND 
vertex. 


Conpling Constants 

= -m|)(GeV 

g'EcNniQ'^ — —m|))(GeV 

Valnes 

12.96 ±3.49 

3.78 ±1.06 


Table 4: Values of the coupling constants g-Zf^NB and gs^ND- 

To sum up, in this work, the strong coupling constants among the heavy bottom spin- 
1/2 Sfe baryon, nucleon and B meson as well as the heavy charmed spin-1/2 baryon, 
nucleon and D meson, namely gj:^NB and gj:^ND, have been calculated in the framework of 
the three-point QCD sum rules. The obtained results can be applied in the analysis of the 
related experimental results at LHC. The predictions can also be used in the bottom and 
charmed mesons clouds description of the nucleon that may be applied for the explanation 
of the exotic events observed by different experiments. These results may also serve the 
purpose of analyzing of the results of heavy ion collision experiments like PANDA at FAIR. 
The obtained results may also come in handy in the determinations of the changes in the 
masses, decay constants and other parameters of the B and D mesons in nnclear medinm. 
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